Solution structure, dynamics and tetrahedral assembly of Anti-TRAP, a homo-trimeric triskelion-shaped regulator of tryptophan biosynthesis in Bacillus subtilis

Graphical abstract


Introduction
The trp RNA-binding attenuation protein (TRAP) senses intracellular levels of tryptophan in order to regulate transcription (Babitzke et al., 1994;Babitzke and Yanofsky, 1993;Gollnick, 1994;Gollnick et al., 1990) and translation (Du and Babitzke, 1998;Merino et al., 1995;Sarsero et al., 2000;Yakhnin et al., 2004Yakhnin et al., , 2001;;Yang et al., 1995) of the genes required for tryptophan biosynthesis (Fig. 1).Transcriptional and translational control by TRAP depends on the ability of tryptophanbound TRAP to bind to an RNA site with multiple (G/U)AG repeats separated by two or three nonconserved spacer nucleotides (Babitzke et al., 1994).The transcription of the six tryptophan biosynthetic genes clustered in the trpECDFBA operon (Henner et al., 1985) is controlled through an attenuation mechanism by which TRAP binding is thought to remodel competing secondary structural elements in the 5 ́ leader region of the nascent mRNA: in the absence of tryptophan, RNA binding by TRAP is impaired, and the leader folds into an anti-terminator hairpin that allows transcription of the entire operon.Increased cellular tryptophan binds to TRAP, activating it to bind to its RNA target which overlaps with and prevents formation of the anti-terminator hairpin; instead, a small downstream terminator hairpin forms and transcription of the structural genes is attenuated (Babitzke and Yanofsky, 1993).There is considerable evidence in support of this terminator/antiterminator model, although TRAP nevertheless induces attenuation in the absence of an intact terminator (McAdams and Gollnick, 2014;Potter et al., 2011).
Translational control of trpE (Du and Babitzke, 1998;Schaak et al., 2003), and perhaps the entire trp operon, is also accomplished through TRAP binding-dependent RNA secondary structure (Yakhnin et al., 2001).When TRAP is inactive and cannot bind to the read-through transcript of the trp operon, the Shine-Dalgarno (SD) sequence of the trpE gene is unstructured and allows ribosome binding.When Trpactivated TRAP binds to the untranslated region of the transcript, the trpE SD sequence becomes sequestered in a hairpin that forms downstream of the TRAP binding site, thereby preventing ribosome binding and translation initiation (Du and Babitzke, 1998).TRAP binding sites are also present in the transcripts of the trpG, trpP, and ycbK genes, but in the translation initiation region, where TRAP binding directly impedes ribosome binding and thereby regulates translation (Babitzke et al., 1994;Du et al., 1997;Lee et al., 1996;Sarsero et al., 2000).
In many bacteria, the extent of charging of tRNA Trp is sensed as a regulatory signal in addition to the levels of free tryptophan (Gutiérrez-Preciado et al., 2009, 2007).In studies with a temperature-sensitive mutant of the tryptophanyl-tRNA synthetase in B. subtilis, the accumulation of uncharged tRNA Trp leads to over-expression of the trpECDFBA operon (Lee et al., 1996).The operon responsible for this effect was shown to be regulated via the T-box transcription anti-termination mechanism, where uncharged tRNA Trp specifically pairs with the leader RNA causing anti-termination and transcriptional read-through (Sarsero et al., 2000).Deletions within the operon identified yczA/ rtpA as the responsible gene, and the protein product, anti-TRAP (AT), was shown to inhibit TRAP activity (Sarsero et al., 2000;Sharma and Gollnick, 2014;Valbuzzi and Yanofsky, 2001).Expression of AT was also shown to be regulated at the level of translation by tRNA-dependent ribosome stalling (Chen andYanofsky, 2004, 2003).
Symmetry mismatch between the oligomeric states of AT and TRAP was evident from early chemical cross-linking, analytical ultracentrifugation, native mass spectrometry and NMR studies (Akashi et al., 2009;Sachleben et al., 2010;Snyder et al., 2004;Valbuzzi andYanofsky, 2002, 2001).Diffracting crystals of Bsu AT revealed a dodecamer (AT 12 ) arranged as a tetrahedral tetramer of trimers (AT 3 ) 4 (Fig. 2) (Shevtsov et al., 2010(Shevtsov et al., , 2005(Shevtsov et al., , 2004)).Each AT protomer forms a short helix (residues 5-8) at the N-terminus, β-strands from residues 9-11 and 34-36 form a short anti-parallel β-sheet, a β-hairpin is formed by strands from residues 20-21 and 24-25, and an α-helix is present at the C-terminus from residues 38-50.Trimers are arranged via a three-helix bundle involving the C-terminal helix, with the zinc binding domains extended outward from the C 3 helical axis.The AT dodecamer is formed through complementary inter-trimer ion pair interactions involving an N-terminal amine and a sidechain carboxylate from each trimer, as well as hydrophobic surface burial of the zinc-binding regions (Fig. 2) (Shevtsov et al., 2010(Shevtsov et al., , 2005(Shevtsov et al., , 2004)).Because the stabilizing ion pair requires a protonated N-terminus, which has a pKa near pH 7 (Sereda et al., 1993), assembly of the dodecamer is disfavored above pH ~ 8 (Sachleben et al., 2010).The B. licheniformis variant of AT (which shares a 76 % sequence identity with B. subtilis AT) crystallized at low pH adopts the same dodecameric structure as the Bsu protein, but at high pH the trimers assemble in an inverted orientation, with their N-termini projecting into solvent (Shevtsov et al., 2010).
AT binds asymmetrically to Trp-loaded TRAP.Although Bsu TRAP is predominantly undecameric, both undecamers and dodecamers are observed in mass spectra (Babitzke et al., 1994;Bayfield et al., 2012;Chen et al., 2011, p. 11;McCammon et al., 2004), and crystallographic studies of Bsu AT bound to dodecameric TRAP rings revealed AT 3 docked on the RNA binding surface of TRAP 12 , via protein-protein surfaces of the zinc binding and helical domains that are not accessible in AT 12 (Watanabe et al., 2009).Mass spectrometric, solution small angle scattering, and analytical ultracentrifugation experiments (Ihms et al., 2014) are also consistent with a model in which only AT 3 can compete with RNA for the TRAP-RNA binding site.These findings suggest that the AT 3 ↔ AT 12 equilibrium may represent an additional regulatory step.
We found experimentally that preventing formation of the intertrimer Met1-Asp7 ion pairs favors AT 3 , even at high protein concentrations (Sachleben et al., 2010).Addition of excess zinc during protein recombinant production of AT in E. coli results in significant accumulation of formylmethionine AT (fAT), presumably due to poisoning of peptide deformylase (Jain et al., 2005;Ragusa et al., 1998); the resulting immature fAT can be chromatographically separated from mature AT under denaturing conditions (Sachleben et al., 2010).Other modification to the N-terminus, such as inclusion of a His-6 tag, similarly prevents AT 12 formation (McElroy, 2005).This provides the means to experimentally produce homogenous samples of either AT 3 or AT 12 , to allow characterization of the two oligomers in the absence of confounding effects from their equilibrium.
We have leveraged the ability to selectively produce AT trimers and dodecamers to determine the solution structure of Bsu AT 3 and to examine the pH-and concentration-dependent equilibrium between to AT 3 and AT 12 using small angle x-ray scattering (SAXS), sedimentation velocity analytical ultracentrifugation (AUC), and native mass spectrometry (nMS).In addition, backbone amide 15 N NMR relaxation measurements (Barbato et al., 1992;Farrow et al., 1994;Kay et al., 1989) allowed us to study the dynamics of both the trimeric and dodecameric AT forms.These results provide new insights into the  structural basis of Anti-TRAP's homo-oligomeric assembly, and how this process may modulate AT's function in its regulatory context.

Anti-TRAP (AT) adopts trimeric and dodecameric species in solution NMR and native MS reveal two states in slow equilibrium
At pH values near 7, purified AT samples at 10 -3 M concentrations yield 1 H-15 N correlation spectra that exhibit doubling of many amide resonances; such doubling arises in cases of chemical heterogeneity or slow conformational exchange between two states.After chromatographically separating mature AT from immature formylated AT (fAT) under denaturing conditions (Sachleben et al., 2010), samples of mature AT produced the same signal doubling (Fig. 3), indicating that it is not due to chemical heterogeneity.Nano electrospray ionization mass spectrometry under soft ionization conditions, at protein concentrations one hundred times lower than NMR (10 -5 M), yielded two sets of ions, one corresponding to the mass of AT 3 , and the other, AT 12 (Fig. 3).Adjusting the solution pH to 8 produced a higher proportion of AT 3 ions, while at pH 6 similar intensities were observed for both sets of ions, consistent with prior findings from pH-dependent NMR spectra.In addition to affecting the oligomeric state, changing the pH slightly altered the preferred charge state distribution: at pH 8, the 7+ and 14+ charge states were favored for AT 3 and AT 12 , respectively; at pH 6, the 8+ and 15+ charge states dominate the spectrum.Although AT 3 and AT 12 populations varied with pH, no other oligomeric species (e.g., AT 1 , AT 6 ) were observed in the nMS data over a wide range of conditions.

Sedimentation velocity experiments provide evidence for a lowly populated hexamer intermediate
Guided by prior sedimentation velocity experiments at a fixed pH that revealed both AT 3 and AT 12 in solution (Snyder et al., 2004), we performed analytical ultracentrifugation sedimentation velocity experiments at pH 7, 7.25, 7.5, and 8 and a range of AT concentrations (Fig. 4).For systems with components undergoing rapid exchange on the timescale of the sedimentation experiment, the modeling of boundary phases can provide the sedimentation coefficients for the components as well as the stepwise binding rates (Brautigam, 2011).Sedimentation velocity data at each pH were globally fit across three concentrations in the program SEDPHAT (Zhao et al., 2015), initially to a simple two-state model comprising the two states observed in NMR and MS spectra, The resulting fits were poor, especially at alkaline conditions, and were better fit with a three-state exchange model, Lamm modeling of the sedimentation data predicts species with an intermediate sedimentation coefficient, which we attribute to a small population of hexamer AT 6 , despite not observing this species by nMS or NMR (Fig. 3).Fits to experimental interference fringe shifts using this model permitted optimization of two sets of apparent rate constants describing dimerization of AT 3 trimers into AT 6 , and a subsequent dimerization of two AT 6 into AT 12 .From the ratio of the fitted rate constants, we obtained consistent equilibrium dissociation constants for the first step, ranging from 20-90 µM over the experimental pH range.For the second step, the association rate constant dropped over three orders of magnitude over the same pH range, resulting in equilibrium constants from 0.05-50 µM.Thus, at low pH the second step is highly favored compared to the first step, resulting in very low concentrations of the hexameric intermediate AT 6 .This is generally consistent with the relatively low populations of intermediate species in the fitted c(s) distributions (Fig. 4).The excellent agreement between experimental and modeled data suggests that monomeric AT is not present at a significant population.

SAXS profiles of AT can be reconstructed from those of AT 3 and AT 12 oligomers
We investigated the pH-dependent structure of AT using small angle x-ray scattering (SAXS).For mixtures, the resulting scattering profiles are the weighted sum of scattering profiles of the individual components (Konarev et al., 2003).Experimental scattering curves were deconvoluted using form factors for trimeric and dodecameric AT (Fig. 5).The predicted scattering profiles for AT 3 and AT 12 were calculated with CRYSOL (Svergun et al., 1995) from the trimer and dodecamer assemblies in the crystal structure (2BX9) (Shevtsov et al., 2005).Deconvolution of the experimental data from the polydisperse AT solutions was performed using OLIGOMER (Konarev et al., 2003).The resulting fits using these two components describe the experimental data well, with AT 12 dominating below pH 7.5.We also tested for the presence of an inverted dodecameric structure such as that observed in crystallographic studies of B. lichenformis AT (PDB 3ld0) (Shevtsov et al., 2010) by including its form factors in fitting the experimental data; this component was found to have zero weight in all of the tested experimental conditions.These data are consistent with a slow exchange between AT 3 and AT 12 under near-physiological conditions, and a negligible population of intermediates.Fig. 2. Oligomeric assembly of Bsu Anti-TRAP (AT).The 53-residue AT protomer (AT 1 ) features a short N-terminal helix leading into a "zinc ribbon" motif and ending in a C-terminal helix.Assembly of AT 3 occurs via a helical bundle that buries ~ 4,200 Å 2 .In solution AT 3 is in an equilibrium with dodecamers (AT 12 ) featuring a tetrahedral arrangement of AT 3 stabilized by a network of ion pairs between the amino terminus of Met1 and the sidechain carboxylate of Asp7 of adjacent trimers.A transient AT 6 intermediate is hypothesized to arise from dimerization of AT 3 via formation of two sets of Met1-Asp7 ion pairs and burial of ~ 1,100 Å 2 ; dimerization of two AT 6 assemblies to form AT 12 is favored by formation of ten additional sets of ion pairs and burial of another ~ 4,600 Å 2 .

NMR studies of fAT 3 Formylated AT is trimeric with high quality NMR spectra
We used immature fAT for detailed investigation of the structure and dynamics of the trimeric TRAP-binding oligomeric form.Since an ion pair between the protonated amine of Met1 and the sidechain carboxylate of Asp7 is proposed to stabilize the dodecameric structure (Sachleben et al., 2010;Shevtsov et al., 2005), the formyl group prevents formation of the ion pair and fAT trimers do not oligomerize to form detectable dodecamers, even at millimolar concentrations at low pH.Resonance assignments and NMR structural data were recorded at 55 • C, to match the conditions used for backbone resonance assignments of Bst TRAP (Kleckner et al., 2012;McElroy et al., 2002;Sachleben et al., 2010).
NMR spectra of fAT 3 are of exceptionally high quality (e.g., Fig. 6), enabling nearly complete backbone and side chain resonance assignments.Bsu AT comprises 53 amino acids, and zinc-bound fAT 3 has a nominal molecular mass of 17.1 kDa.Three-fold symmetry results in spectra reflecting the number of signals of a single chain, while fast tumbling result in narrow NMR resonances and high-quality heteronuclear correlation spectra.Backbone resonances were assigned from standard triple resonance spectra, while sidechain assignments were obtained from a combination of 2D and 3D double-and triple-resonance spectra.This allowed for assignment of every proton in the protein except for the side chain amine and hydroxyl protons.
The structure of fAT 3 is unperturbed by pH SAXS, nMS and NMR data demonstrate that the trimer-dodecamer equilibrium is strongly affected by pH (Fig. 3, Fig. 4, Fig. 5).This behavior was attributed to a near-neutral pK a for the amine at the Nterminus of the protein (Sereda et al., 1993), and its ion pair with the carboxylate side chain of Asp7 (Sachleben et al., 2010;Shevtsov et al., 2005).There is one other functional group in AT with a predicted pK a near neutrality: His50, which is conserved among AT variants (Shevtsov et al., 2010;Valbuzzi and Yanofsky, 2001).In crystals the side chain imidazole Nε2 of His50, located at the end of the C-terminal helix, is observed to form an intra-chain ion pair with the carboxylate of Glu9, located in the zinc-ribbon motif (Fig. 2).Changes in the protonation state of the imidazole side chain could conceivably alter the geometry between the zinc-binding and helical domains and thus alter trimer packing in the dodecamer, driving changes in oligomeric structure (e.g., such behavior has been engineered into homo-oligomeric assemblies (Boyken et al., 2019)).To examine whether His50 or another titratable group in the protein might induce structural changes that contribute to the pH-dependent oligomerization behavior, we recorded 1 H-15 N correlation NMR spectra of fAT over a range of pH values (Fig. 7).Since fAT has a blocked N-terminus and does not form stable dodecamers, this allowed distinguishing spectral effects of ionization from those resulting from oligomerization.Correlation spectra of fAT recorded over a pH range that would shift the predominant oligomeric state from AT 3 to AT 12 in the mature protein result in only minor shifts in the spectra, particularly the amides of Lys49 and Leu51, likely reflecting their proximity to the His50 indole.This is consistent with a dominant role of a near neutral pKa for the N-terminal amide group in stabilizing the dodecamer.

Solution structure determination
The solution structure of the TRAP-binding oligomeric state of Bsu AT (AT 3 ) might be inferred from the 2.8 Å crystal structure of the dodecamer (2BX9) (Shevtsov et al., 2005), the 2.1 Å structure of the Bli dodecamer (3LCZ) (Shevtsov et al., 2010), and 3.2 Å structure of the Bsu AT bound to Bst TRAP (2ZP8) (Watanabe et al., 2009), but de novo structure determination was justified by the potential for perturbations from oligomerization and crystal packing.To that end, we determined the solution structure of fAT 3 using distance restraints from NOEs, orientation restraints from residual dipolar couplings, dihedral angle restraints from three-bond scalar couplings, three-fold symmetry restraints, and overall shape restraints from SAXS scattering curves.
Because of its homo-oligomeric structure and diverse source of restraints, structure determination proceeded via several steps.First, three-dimensional models of the zinc-bound AT protomer were computed from torsion angle restraints and NOEs that could be unambiguously assigned as being intramolecular.Next, we determined the solution structure of the trimer by replicating the chain and assigning as ambiguous all NOEs for which conservative distance filters prevented their assignment as being intramolecular.Most ambiguous NOEs could be assigned via iterative structure refinement steps, though some NOEs remained ambiguous throughout the refinement process (Table 1).
Symmetry, residual dipolar coupling (RDC) and SAXS restraints significantly increased the precision of the NMR structure (Fig. 8).Ensembles computed using only NMR-derived distance and torsion angle restraints yielded poor convergence, with an overall backbone RMSD of 1.5 Å, and per-protomer backbone RMSD of 1.3 Å. RDCs and SAXS data improved the overall precision, but not as much as inclusion of threefold symmetry restraints, which increased the overall precision by almost an Ångstrom (from 1.5 Å to 0.6 Å).

Comparison to crystallographic and AlphaFold2 models
The final NMR ensemble was well defined by the restraints, with a backbone RMSD of 0.3 Å for C α atoms (Fig. 8, Fig. 9, Table 1), and reveals a high degree of similarity with the trimers in the B. subtilis dodecamer (RMSD ~ 0.6 Å; Fig. 9).This is similar to that observed between protomers in the 2.8 Å crystal structure of AT 12 (~0.6Å for backbone atoms in 2BX9), which was determined with 12-fold noncrystallographic symmetry (Shevtsov et al., 2005(Shevtsov et al., , 2004)), and within the variance between those structures and that predicted by AlphaFold2 (Jumper et al., 2021).We used the longitudinal and transverse NMR relaxation times (T 1 , T 2 ; the inverse of the relaxation rates R 1 , R 2 ) of backbone amide 15 N nuclei to validate the solution structure (Tjandra et al., 1997).In the absence of internal motions (see below) these values are dominated by reorientation of the amide nitrogen-proton bond vector and reflect the axially symmetric diffusion of the trimer.The experimentally measured T 1 /T 2 ratios for AT 3 are in good agreement with those predicted from members of the NMR ensemble, and from the trimer in the AT 12 crystal structure (Fig. 9) (Bernadó et al., 2002;García de la Torre et al., 2000).The largest deviations between the measured and predicted values arise at the N-and C-termini and the zinc binding loops, suggesting both subtle differences between the structure in crystals and in solution, as well as fast and slow conformational averaging in solution that is not possible in the crystal lattice.

AT dynamics
Backbone amide 15 N R 1 , R 2 , and 1 H-15 N NOE relaxation rates that are sensitive to reorientation of the amide 1 H-15 N bond vector (Kay et al., 1989) were measured on both the 68 kDa dodecameric AT 12 and 17 kDa trimeric AT 3 using conventional two-dimensional heteronuclear correlation spectra (Fig. 10) (Barbato et al., 1992;Farrow et al., 1994).Differences in the R 1 and R 2 relaxation rates for the two oligomeric states reflected the large difference in their rotational diffusion tensors; overall correlation times τ m computed from the 10 % trimmed mean R 2 /R 1 ratios were 8 and 25 ns at 55 • C, respectively (Lee et al., 1997).Hydrodynamic calculations based on the crystal structure of AT 12 , and the R 2 /R 1 ratios were also consistent with an isotropic diffusion tensor and an overall τ m of 25 ns.However, R 2 /R 1 ratios for AT 3 were consistent with an axially symmetric diffusion model with an anisotropic diffusion tensor, D iso = 2.00 ± 0.01 x 10 7 s − 1 , D ‖/⊥ = 0.82 ± 0.04 (τ m = 8.3 ns).
Large differences in the rotational diffusion tensors of AT 12 and AT 3 are accompanied with differences in internal motions.Amide 15 N relaxation rates were analyzed using the extended model-free formalism (Clore et al., 1990;Lipari andSzabo, 1982a, 1982b) to compare the amplitude internal dynamics of AT 3 and AT 12 on the ps-ns timescale, as represented by the square of the generalized order parameter, S 2 (Fig. 11).While both proteins exhibited generally high-order parameters, AT 3 and AT 12 were found to exhibit local differences in fast internal motions.High S 2 values for AT 3 imply rigidity of the C-terminal helical bundle, and larger amplitude fast motions in the zinc binding arm domains and the N-terminal helix.By contrast, the zinc binding domains in AT 12 are relatively rigid on the ps-ns time scales reported by S 2 , likely a reflection of their packing onto adjacent zinc binding domains in the dodecamer (Fig. 2).The two oligomeric forms of AT also exhibit different patterns of the phenomenological exchange broadening term R ex , attributable to dynamics on slower µs-ms time scales (Clore et al., 1990;Mandel et al., 1995).For AT 3 , the residues with largest R ex terms are located at the N-terminus of the protomers, with smaller R ex terms at the zinc binding domains.For AT 12 , almost the inverse is observed, with the largest exchange terms needed to fit data for the helical bundle.

Discussion
Tryptophan biosynthesis is highly regulated in bacteria (Gutiérrez-Preciado et al., 2007;Merino et al., 2008;Yanofsky, 2007).Within bacteria, TRAP-mediated transcription attenuation and translation repression via RNA binding represents an unusual regulatory mechanism.In the subset of organisms that also encode Anti-TRAP, its expression represents yet an additional regulatory layer by binding to and reversing the activity of an inhibitory protein.In this work we examined the equilibrium between active and inactive oligomeric states of Bsu Anti-TRAP using nMS, SAXS, AUC and NMR spectroscopy.We determined the solution structure of the active AT 3 form of the protein and observed differences in its internal dynamics relative to that of the inactive form AT 12 .These studies thus help to inform our understanding of how AT inhibits RNA binding by Trp-bound TRAP, and of how the AT 3 ↔ AT 12 equilibrium may provide yet another regulatory layer.
In our studies of Bsu AT, serendipity played a significant role.The fortuitous requirement of a mature amino-methionine terminus for stabilizing the AT dodecamer coincided with experimental zinc supplementation of minimal media for recombinant protein expression, which likely prevented maturation of a sizeable fraction of the protein.
The second residue in Bsu AT is a valine, which is often a poor substrate of methionine aminopeptidase (Frottin et al., 2006;Hirel et al., 1989), thus ensuring retention of Met1 in the mature protein.The application of NMR and MS enabled detection of the formylated N-terminus, while denaturing reversed-phase chromatography allowed separation of immature fAT from mature AT, which was not possible with conventional native purification.While the AT 3 ↔ AT 12 equilibrium is evident in samples of mature AT, the ability to prepare complexes of pure fAT provided an unambiguous means of establishing the role of the N-terminal amine in the equilibrium and enabling independent structural characterization of the TRAP-binding form of the protein (AT 3 ).
We found that Bsu fAT 3 (and AT 12 ) yields NMR spectra of exceptional quality, enabling complete resonance assignments and high-resolution structure determination by application of NMR, symmetry and SAXS restraints.Symmetry restraints were found to strongly increase the precision of the structural ensemble, while SAXS and RDC restraints provided modest improvements.We found the solution structure of AT 3 to be indistinguishable from that in the 12-mer crystal, or from the AF2 predicted model, indicating little conformational adjustment upon formation of the dodecamer.
While the average structure of the AT trimer in isolation and in the dodecamer are highly similar, they differ in their dynamics.Consistent with the overall size and shape of the oligomers, backbone 15 N relaxation measurements revealed an isotropic rotational diffusion tensor for AT 12 and axially symmetric diffusion for the smaller AT 3 .Model-free analysis of internal dynamics revealed flexible zinc ribbon motifs in AT 3 , suggesting a role for the dynamics of these motifs in engaging Trpbound TRAP rings.AT 12 has comparatively more rigid zinc-binding domains, likely reflecting packing constraints from the other trimers in the tetrahedral arrangement.Modest R ex terms in the N-terminal and zinc-binding domains of AT 3 may reflect exchange dynamics on the µs-  We conclude that the Bsu AT 3 -AT 12 equilibrium is dependent on the near-physiological pKa of the Met1 amine, with AT 12 favored by a protonated N-terminus and AT 3 favored when unprotonated.Prior NMR observations (Sachleben et al., 2010) are corroborated here by nMS, SAXS, and AUC experiments as a function of pH.NMR spectra of fAT, over a pH range that alters the trimer-dodecamer equilibrium for mature AT, only show minor shift perturbations adjacent to His50, supporting the dominant role of the protonated N-terminus in stabilizing the dodecamer (Fig. 7).Native MS spectra (at lower protein concentration) reveal AT 3 as the major species at pH 8, while at pH 6 nearly equal intensity is observed for ions arising from AT 3 and AT 12 species (Fig. 3); this provides some context to a prior nMS study which found exclusively AT 12 when sprayed from pH-unadjusted ammonium acetate (Akashi et al., 2009;Watanabe et al., 2009), which can approach a value of 4.75 in the ESI plume (Konermann, 2017).SAXS data for mature AT recorded as a function of pH and protein concentration could be well described by considering varying ratios of AT 3 and AT 12 (Fig. 5).Experimental curves Fig. 6.Symmetry and exceptional NMR spectra of fAT 3 allowed nearly complete resonance assignments.Aliphatic region of a two-dimensional constant-time 1 H-13 C correlation spectrum of fAT 3 .In this spectrum 13 C resonances coupled to even (red) and odd (black) numbers of aliphatic carbons have opposite sign.Aliphatic carbon resonance assignments as indicated.could be reproduced form a linear combination of the form factors computed from the crystallographic dodecamer and a trimer extracted from the dodecamer.This finding is consistent with the failure to observe other assembly intermediates under equilibrium conditions.Sedimentation velocity (AUC) experiments further corroborate the pHdependent equilibrium between AT 3 and AT 12 and provide evidence for a short-lived AT 6 intermediate.If we consider AT 3 a "monomer" that can tetramerize, this behavior is consistent with biophysical characterization of other "tetramers", which predominantly follow a monomer--dimer-tetramer oligomerization path (Powers and Powers, 2003).Small surface burial and minimal polar contacts in AT 6 , result in transient/weak dimerization; large surface burial and many polar contacts favor AT 12 upon collision of two AT 6 "dimers of trimers".

Conclusion
The mechanism by which Anti-TRAP (AT) promotes Trp production by inhibiting RNA binding by tryptophan-activated TRAP remains poorly understood (Ihms et al., 2014;Snyder et al., 2004;Valbuzzi and Yanofsky, 2001;Watanabe et al., 2009).The capacity of AT to form stable dodecameric and trimeric structures, wherein only the trimeric form is competent for binding and inhibiting TRAP is suggestive of a regulatory role for that equilibrium.While the complementary biophysical studies reported here provide qualitatively consistent insights into the properties of AT and its mechanism of oligomerization and TRAP inhibition, further work would be required to establish quantitative consistency.Some sources of uncertainty include the effect of molecular tumbling on observed NMR signal intensities, of ionizable buffers and of protein charge on gas phase ionization and gas phase signal detection, and relatively low structural resolution of SAXS and sedimentation velocity methods.Those caveats aside, it is clear that in vitro AT exhibits a slow equilibrium between active AT 3 and inactive AT 12 states.
This slow association/dissociation behavior of AT could serve a regulatory function by kinetically modulating the response to changes in Trp concentration (Frieden, 2019).In this model (Fig. 1), TRAP is slow to release Trp upon decrease in its cellular concentration and remains active for RNA binding and inhibiting trp expression.Meanwhile, decreasing cellular Trp results in accumulation of uncharged tRNA Trp and T-box mediated activation of AT expression.Expressed AT will trimerize and rapidly bind to RNA-free Trp-activated TRAP rings, preventing them from binding new trp RNA leaders.AT 3 that do not find TRAP rings could instead dodecamerize, provided the cellular pH is near the pKa, resulting in slowly equilibrating populations of active AT 3 and inactive AT 12 .This scenario would thus generate a pool of fast-acting AT 3 and a reservoir of slowly activated AT 12 .The resulting steadystate equilibrium between AT 3 , AT 12 , AT 3 -Trp-TRAP 11 , Trp-TRAP 11 and Trp-TRAP 11 -RNA would achieve a steady state of Trp homeostasis.An alternative intriguing possibility is that the pH-dependent AT 3 -AT 12 equilibrium serves to regulate gene expression in response to changes in intracellular pH.AT is not known to be associated with bacterial response to changes in pH, but sporulation by B. subtillis is a consequence of changes in gene expression patterns and is associated with decreases in cell volume and pH that would likely affect the AT 3 -AT 12 equilibrium.The notion that pH-driven changes in AT oligomeric state could serve as a mechanism for sensing and responding to such environmental changes is an intriguing possibility that merits further exploration (Errington, 1993;Magill et al., 1994).
The properties of the AT 3 "triskelion", with its polydentate configuration of flexible zinc-binding arms, seem uniquely suited for capturing one or multiple TRAP rings, and thereby preventing RNA wrapping.The use of the amino terminus to mediate oligomerization represents a uniquely tuned equilibrium since the pKa and association constants are coupled, such that high protein concentration would shift the pKa to favor the protonated state, and thus the oligomer.It remains unclear why some species of Bacilli evolved to produce AT as an extra layer of regulation (e.g., Bsu, Bli), while others did not (e.g., Bst, Bha).This study highlights the lengths to which organisms have evolved sophisticated mechanisms for metabolic control, and the power of complementary biophysical measurements to describe those properties that enable regulated function.

Growth and purification of B. subtilis Anti-TRAP
AT was recombinantly expressed from E. coli from an engineered codon-optimized expression vector (Sachleben et al., 2010).B. subtilis Anti-TRAP is a 53 amino acid polypeptide, with the sequence: MVIATDDLEV ACPKCERAGE IEGTPCPACS GKGVILTAQG YTLLDFIQKH LNK.It has a molecular weight of 5649.57Da and its predicted extinction coefficient at 280 nm is 1490 M − 1 cm − 1 in its reduced form (Expasy ProtParam).Uniformly N-terminal formylmethionine Anti-TRAP (ƒAT) and deformylated AT were purified as previously described (Sachleben et al., 2010).Briefly, transformed BL21(DE3) cells grown and induced in modified M9 minimal media were lysed and cellular debris pelleted by centrifugation.Isotope labeling was achieved by replacing unlabeled ammonium chloride with 15 N ammonium chloride for 15 N labeling, by substitution of 13 C glucose for 13 C/ 15 N double labeling.Heating of the cell lysate supernatant for 15 min at 70 • C and subsequent centrifugation removed most native E. coli proteins.The resulting supernatant was loaded onto Q-sepharose beads (GE healthcare) and eluted via a NaCl gradient.Eluted Anti-TRAP was denatured by addition of 6 M guanidinium and subjected to reverse-phase HPLC purification.AT (i.e., deformyl) and ƒAT elute at sufficiently resolved retention times from a C4 column to be separated and lyophilized, providing > 99 % desired species.Subsequent refolding results in intact ƒAT trimers for ƒAT, and AT 3 + AT 12 in equilibrium for (deformylated) AT.

Mass spectrometry
The theoretical mass of AT (from ExPASy ProtParam): Monomer (without Zn) = 5,649.57Da; monomer with Zn (65.38 Da) = 5714.95Da; trimer with three Zn = 17,144.85Da; dodecamer with twelve Zn = 68,579.40Da.AT samples were buffer exchanged twice into 100 mM ammonium acetate at 3 different pH's (pH 6.8, unadjusted; pH 8, adjusted using ammonium hydroxide; pH 6, adjusted using acetic acid) using size-exclusion spin columns (MicroBioSpin-6; Bio-Rad); the final concentration of AT was 10 μM.Borosilicate glass capillary emitters were pulled in-house with a Sutter Flaming/Brown Micropipette puller P97; samples were loaded and sprayed using nanoESI.Mass spectra were recorded on a Waters Synapt G2S mass spectrometer modified in-house with a SID device located between the trap cell and the drift cell.Data were recorded with capillary voltage, 0.7-1.0kV; cone voltage, 20 V; extraction cone, 5-8 V; helium gas flow, 120 mL/min; trap gas flow, 2 mL/min; ion mobility nitrogen gas flow, 60 mL/min; source temperature 20 • C; trap wave velocity, 160 m/s; trap wave height, 4 V; ion mobility wave velocity, 300 m/s; ion mobility wave height, 20 V; transfer wave Fig. 9. Validation of the AT 3 solution structure.(A) Solution structure of AT 3 is congruent with its structure in AT 12 crystals and that predicted by AlphaFold2.Left, green, NMR ensemble of 20 conformers for chain A; grey, superposition of the twelve chains in the 2BX9 crystal structure; pink, AlphaFold2 prediction.Right, superposition of the AT 3 NMR derived ensemble with chains A-C (blue, orange, green) in the 2BX9 crystal structure.(B) The ratios of the NMR T 1 and T 2 relaxation times for backbone amide 15 N provide an independent validation of the NMR-determined structure.Blue, experimental data; green, predicted ratios based on hydrodynamic modeling of the top five models of the NMR ensemble; black, predicted ratios using chains A-C from the crystal structure.Secondary structure indicated schematically; "*" mark the zinc ligands.

NMR assignments and structural restraints of AT 3
fAT spectra exhibit signals exclusively attributable to trimeric AT.NMR spectra of fAT were recorded on a 3 mM sample in 20 mM d 11 -Tris pH 7.5, 20 mM β-mercaptoethanol, 1 mM ZnCl 2 and 0.02% (wt/vol) NaN 3 The NMR spectra of fAT3 also exhibit the number of peaks expected for the monomer, consistent with its molecular symmetry.Assignments of uniformly, doubly labeled ( 13 C, 15 N) trimeric AT were obtained using HNCO, HNCACB, CCONH-TOCSY, HCCONH-TOCSY, HCCH-COSY, HCCH-TOCSY, 15 N-edited TOCSY, constant-time 13 C-edited HSQC, 13 C-edited HSQC, constant-time 13 C-edited HSQC of the aromatic region, 13 C-edited HSQC of the aromatic region, and a 1 H-1 H plane and 1 H-13 C plane of the HCCH-COSY of the aromatic region (Cavanagh et al., 2006).
Torsion angle restraints for the backbone φ angles were obtained from the analysis of 3D HNHA spectra (Bax et al., 1994) and additional backbone torsion angle φ, ψ restraints were obtained from chemical shifts using TALOS (Cornilescu et al., 1999).Sidechain χ 1 torsion angle restraints were obtained from the HNHB and HACAHB-COSY experiments (Bax et al., 1994).The X-Pro amide bonds for AT 3 were determined to be in the trans conformation from analysis of the Cβ and Cγ chemical shifts (Schubert et al., 2002); we did not assign Cγ resonances for AT 12 but all amides are modeled as trans in the crystal structure (Shevtsov et al., 2005).One bond 1 D HN-N and 1 D N-C residual dipolar coupling (RDC) restraints were obtained from IPAP spectra (Ottiger et al., 1998;Permi, 2003), of 13 C/ 15 N-labeled AT aligned with 20 mg/ml Pf1 phage (ASLA Biotech).NOE restraints were obtained from a 13 C/ 15 N-edited time-shared NOESY (Sattler et al., 1995) in H 2 O and a 13 C-edited NOESY (Cavanagh et al., 2006) in 2 H 2 O.

NOE assignments and initial AT 3 structure determination
NOE assignments of the symmetrical homo-trimer were made using the symmetry ambiguous distance restraint (symmetry ADR) method starting with a well-defined monomer (Kovacs et al., 2002;O'Donoghue et al., 2000).The structure of the protomers were first obtained by assigning NOEs presumed to be intramolecular, (i.e., NOEs between zinc-coordinating coordinating cysteines, secondary structural elements inferred from secondary structure predictions based on sequence and chemical shift data), and using these as input for restrained torsionangle refinement with CNS (Brunger, 2007).The set of 20 lowest energy structures (of 35 total) was then used with a 10 Å cutoff to filter through the remaining NOE data.At this stage, NOEs with less than six ambiguous assignments were used for another round of structure refinement and the cycle was iterated while decreasing the distance cutoff by 2 Å until the cutoff was 6 Å and the monomer structure converged to an RMSD of 2.08 Å for backbone atoms.The trimer structure was then calculated by ambiguously assigning any NOEs that could not be satisfied intramolecularly (i.e., they were defined as being between protomers A, B and C) and a structure calculation was performed utilizing non-crystallographic symmetry restraints in addition to the distance, orientation, and torsion angle restraints.This set of structures was then used with a 6 Å filter to assign the remaining NOEs and sort out the interprotomer contacts.

SAXS and RDC refinement of AT 3
To obtain a scattering profile for the trimer, pure ƒAT was purified, refolded and dialyzed against 50 mM sodium-phosphate, pH 8, 100 mM NaCl, 200 µM ZnCl 2 and 1 mM TCEP.SAXS profiles at 1, 2, 3.3, 5, and 10 mg/mL ƒAT (monomer) were acquired at the SIBYLS beamline at the Advanced Light Source in Berkeley, California.Data were obtained at a Q spacing range of 0.01 to 0.32 Å-1, and at 0.5, 1, and 6 s exposure times.The scattering contribution from the buffer was removed by subtracting the scattering profile of the dialysis solution.Processing, validation, and analysis of SAXS data was performed with the ATSAS package (Manalastas-Cantos et al., 2021) and custom scripts.For restraints to the AT 3 structure, the 6 s exposure profile of the 3.3 mg/mL sample was used, as a small concentration-dependent amount of dodecamer (<3% at 10 mg/mL) was detectable from deconvolution of the higher concentrations.
Isotropic and anisotropic IPAP spectra for RDC data were collected in 50 mM sodium-phosphate pH 6.6, 2 mM DTT at 298 K in the absence and presence of 20 mg/mL Pf1 phage as the alignment medium.The lowest energy initial trimer structure (from above) was used as the starting point for simultaneous refinement against NOE, dihedral, RDC, and SAXS restraints in Xplor-NIH (Schwieters et al., 2018).Agreement between measurements and values calculated from the structure was assessed using the DC server (https://spin.niddk.nih.gov/bax/nmrserver/dc/).Following refinement, the AT 3 structure ensemble was validated using protein structure validation server (PSVS) (Bhattacharya et al., 2007) and Procheck-NMR (Laskowski et al., 1996).Final refinement and validation statistics are provided in Table 1.

Rotational diffusion
Overall rotational correlation times, τ m , and rotational diffusion tensors were determined experimentally from R 2 /R 1 ratios (Lee et al., 1997) utilizing only R 2 /R 1 ratios for residues with relaxation rates within one standard deviation of the 10 % trimmed mean value.These calculations for AT 12 showed no statistical significance for fitting the NMR data to an axially symmetric or fully anisotropic diffusion tensor over an isotropic model.Therefore, an isotropic diffusion tensor with a τ m of 25 ns was used for the model-free calculations of the dodecamer.This was supported by calculations of the hydrodynamic parameters by HYDRONMR 7.C (Bernadó et al., 2002;García de la Torre et al., 2000) based on the dodecameric structure 2BX9 that also predict nearly isotropic tumbling for AT 12 ; the crystallographic coordinates were protonated using PDB2PQR (Jurrus et al., 2018).The calculations based on the R 2 /R 1 ratios of AT 3 , however, indicated a statistically significant improvement for the axially symmetric model over the isotropic model.This was similarly supported by hydrodynamic calculations that yielded an anisotropic diffusion tensor with two similar and one unique axis.Therefore, an axially symmetric diffusion tensor from the R 2 /R 1 ratios was used for the model-free calculations of AT 3 .Validation of the resulting structures from the R 2 /R 1 ratios was performed with HYDRONMR 7.C, using a solvent viscosity of 0.005P, T = 328 K, 1.02 Å N-H bond distance, CSA of − 160 ppm and 14.09 T field.Diffusion tensors were computed for the symmetric trimers and T 1 /T 2 ratios were analyzed for the first chain (Fig. 9).

Model-free analysis
Model-free analysis for both AT 3 and AT 12 (Clore et al., 1990) was performed using the program TENSOR v. 2 (Dosset et al., 2000), which was recompiled to allow for input of up to 1000 residues and a CSA of − 172 ppm (Lee and Wand, 1999).Five motional models with up to three parameters describing the internal motions were used: (1) S 2 (2) S 2 and t e (3) S 2 and R ex (4) S 2 , t e , and R ex (5) S 2 , S 2 f , and t e .Standard errors in the model-free parameters were obtained from Monte Carlo simulations using 1000 samples.Model fitting was performed using TENSOR v. 2 via a modified step-up hypothesis testing method (Mandel et al., 1995).

Observation of the pH-dependent AT 3 -AT 12 equilibrium via SAXS
The program CRYSOL (Svergun et al., 1995) was used to generate predicted scattering profiles AT 3 structure, the B. subtilis AT dodecamer retrieved from the RCSB (3BX9), and the B. licheniformis high-pH "inverted" dodecameric isoform (3LD0).Given the structural similarity of the AT trimer component between B. subtilis and B. licheniformis at nominal SAXS resolution (10-20 Å), the B. licheniformis dodecamer coordinates were used without modification.All three profiles were normalized to the square of the input structure's molecular weight.The solvent density and contrast of the hydration shell were determined by fitting the scattering of the AT trimer structure against an empirical scattering profile of 5 mg/mL fAT in CRYSOL.The empirical scattering data from mature AT at 16 different pH and concentration conditions was fit with the program OLIGOMER (Konarev et al., 2003), using the three predicted scattering profiles.An additional systematic term was also used to correct for baseline offset arising from incomplete buffer subtraction/mismatch.

Analytical ultracentrifugation
AT stock solutions were split and dialyzed against sample buffers at pH 7, 7.25, 7.5, and 8, consisting of 50 mM sodium phosphate, 100 mM NaCl, 0.5 mM tryptophan, and 0.02 % (w/v) NaN 3 .Sedimentation samples were subsequently made from the concentrated stock solutions with the corresponding dialysis buffer and allowed to equilibrate overnight.To account for the decreasing self-affinity at increasing pH, samples at pH 7-7.25 were made to 250, 125, and 65 µM AT monomer concentration, and at pH 7.5-8, to 500, 250, and 125 µM.Dilutions were performed on a gravimetric balance accurate to 0.2 mg, with an assumed solution density of 1.0 g/ml.Sedimentation velocity data were obtained at 25 • C using a Beckman Coulter ProteomeLab XL-I ultracentrifuge equipped with an 8-position rotor and double-sector cells with sapphire windows.Matched dialysis buffer was present in the reference sector.Data were recorded using interference optics with a scan interval of one minute and 200 scans per cell at a rotor speed of 50,000 rpm.The position of the meniscus and cell bottom were determined for each cell by fitting a sedimentation coefficient distribution and uniform frictional value f/f 0 in SEDFIT (Schuck, 2013) at a sampling interval of 10 to 20 points per Svedberg.Time and radial invariant noise were removed from each interference scan.The buffer density ρ and viscosity η under these conditions were calculated by SEDNTRP (Philo, 2023) to be 1.01260 g/ mol and 1.0520 cP, respectively.
Determination of kinetic parameters for the AT 3 ↔ AT 6 ↔ AT 12 model in SEDPHAT (Zhao et al., 2015) was performed by globally fitting 200 interference scans from each set of AT concentrations at one pH (600 scans total).In addition to model kinetic parameters, the fringe shift coefficient, apparent sedimentation coefficients for AT 3 , AT 6 , and AT 12 , as well as the molar mass for the trimer were also optimized.pH 8 data were fit first, using the HYDROPRO (García De La Torre et al., 2000) -predicted sedimentation coefficients of 2.1 S and 4.1 S at 25 • C for AT 3 and AT 12 , respectively, as starting sedimentation coefficients.The AT 6 starting sedimentation coefficient was set to 3 S. Initial parameters were 10 5 M − 1 for k 1 , 10 -3 s − 1 for k -1 , 10 6 M − 1 for k 2 , and 10 -4 s − 1 for k -2 , as these provided reasonable predicted interference curves.The apparent molar weight of AT in the experiment was also optimized, using the expected mass of the trimer as a starting value.A reasonable starting point for the molar absorptivity was determined by manual iterative evaluations of the model using different values until a fringe shift of appropriate magnitude was achieved.(Uncertainty analysis in the kinetic parameters was attempted through Monte-Carlo estimation of errors, but after 1000 evaluations no variance in the fitted parameters was detected.).

Interfaces
Solvent accessible surface (SAS) area burial between protomers was computed using PyMOL 2.5.2 (Schrödinger, LLC) using a dot density of 3. Burial was computed from difference between the SAS of components and their complex, using PyMOL get_area.SAS burial for folding of the protomer was estimated from multiple conformations using PROSA (Estrada et al., 2009).Similar interface and surface areas were obtained using PISA (Krissinel and Henrick, 2007).
Coordinates of fAT 3 have been deposited to the RCSB, accession 9CCH, and NMR data, including assignments and restraints obtained 328 K deposited to the BMRB, accession 16492.Raw data, including NMR spectra, SAXS curves and sedimentation velocity data are being deposited to BMRBig, accession (bmrbig100).

Fig. 1 .
Fig. 1.Transcription and translation of the trp operon in B. subtilis is regulated by feedback loops involving the functions of two proteins: the ring-shaped 11-mer TRAP, which senses free Trp, and Anti-TRAP (AT) which is overexpressed in response to accumulation of uncharged tRNA Trp .In the absence of AT, Trp-bound TRAP (blue) binds a series of UAG and GAG triplets in the 5′ leader of the trp operon, resulting in transcriptional and translational repression in part by formation of a terminator or sequestering the Shine-Dalgarno (SD) sequence.When expressed, trimeric AT (AT 3 ) can bind to Trp-activated TRAP, preventing RNA binding and blocking its inhibition of Trp biosynthesis.AT exists in equilibrium between trimeric (red) and dodecameric states (blue, orange, green; the red trimer is behind the other three), while only trimeric AT binds and inhibits TRAP.

Fig. 4 .
Fig.4.Sedimentation velocity experiments over a range of pH values can be fit with a three-state association model AT 3 ↔ (AT 3 ) 2 ↔ AT 12 .Interference fringe shifts (gray dots) are globally fit across three AT concentrations (first three columns).Fits provided by the mechanistic model are in blue, green, and red lines.Only every tenth scan is shown.Residuals are below the data.Fourth column, sedimentation coefficient distributions with blue, red and green corresponding to the concentrations on the left, at pH 7.00 (A), pH 7.25 (B), pH 7.50 (C), pH 7.75 (D), and pH 8.00 (E).C.A. McElroy et al.

Fig. 5 .
Fig. 5. SAXS profiles of AT can be described as mixtures of AT 3 and AT 12 .A, scattering curves for AT 3 and AT 12 , computed from the dodecameric crystal structure (2bx9).B-D, left, experimental scattering intensities (black) for AT recorded over a range of pHs (6, 7, 7.5, 8) and concentrations (1.77, 0.89 and 0.44 mM), fit to a two-state model with varying fractions of AT 3 and AT 12 .Profiles are offset in the y-axis from each other for visibility.Right, fitted volume fractions as a function of pH for each sampled concentration.

Fig. 7 .
Fig. 7. Protonation of the Met1 amine is coupled to oligomerization of AT.Left, 1 H-15 N correlated spectra of AT recorded at 25 • C and pH 6.5, 7.0, 7.5 and 8.0 reveal strong perturbations of several signals, and line narrowing at high pH, which favors trimeric AT.Right, for fAT the N-terminal amine is protected by the formyl group and very minor shift perturbations are observed only for residues flanking His50; this indicates no pH-dependent structural shifts in the AT trimer.Amide assignments transferred from 55 • C via a temperature series.(sc), side chain amide.
-restricting constraints and structure quality factors from protein structure and validation server (PSVS).

Fig. 8 .
Fig. 8. Impact of restraint classes on precision of AT 3 NMR ensemble.Each ensemble was computed using NMR-derived torsion angle restraints, and (A) NOEs, (B) NOEs and three-fold symmetry restraints, (C) NOEs, symmetry, special restraints from SAXS data, (D) NOEs, symmetry restraints and residual dipolar couplings, and (E) all restraint classes.Backbone RMSDs as indicated.

Fig. 10 .
Fig. 10.Amide 15 N relaxation rates for trimeric and dodecameric AT reveal differences in overall correlation and local dynamics.Relaxation rates were measured at 600 MHz, 55 • C. Local correlation times were computed from the R 2 /R 1 ratios.Secondary structures and location of metal ligands are indicated schematically above.

Fig. 11 .
Fig. 11.Fitted Model-Free dynamics parameters suggest a stable helical bundle and flexible zinc-binding domains in AT 3 , but less variable dynamics in AT 12 .Top, fitted model free parameters S 2 and R ex for AT 3 and AT 12 .Middle, fitted S 2 values mapped to the backbone of the AT trimer using a linear color ramp from the mean value (grey) to + (blue) or -(red) two standard deviations from the mean.Bottom, R ex values mapped to the AT 3 structure with a color ramp from 0 to ≥ 4/ s (yellow).

Table 1
NMR structural statistics for AT ensemble (20 structures):